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ABSTRACT 

We report the discovery of a new double image gravitational lens system B1030+074 
which was found during the Jodrell Bank - VLA Astrometric Survey (JVAS). We have 
collected extensive radio data on the system using the VLA, MERLIN, the EVN and 
the VLBA and optical observations using WFPC2 on the HST. The lensed images are 
separated by 1.56 arcseconds and their flux density ratio at centimetric wavelengths is 
approximately 14:1 although the ratio is slightly frequency dependent and the images 
appear to be time variable. The HST pictures show both the lensed images and the 
lensing galaxy close to the weaker image. The lensing galaxy has substructure which 
could be a spiral arm or an interacting galaxy. 

Key words: gravitation - galaxies: individual: B1030+074 - gravitational lensing. 
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1 INTRODUCTION 



The Jodrell-Bank VLA Astrometric Survey (JVAS) is a sur- 
vey of flat-spectrum radio sources one of whose purposes is 
to search for gravitational lens systems. JVAS contains the 
~2400 strongest flat spectrum sources in the northern sky 
(Patnaik et al. 1992a; Patnaik 1993a; Browne et al. 1998, 
Wilkinson et al. 1998). All sources were observed with the 
VLA at 8.4 GHz giving a resolution of 200 mas. Those that 
were found to possess multiple components or very com- 
plex structure have been followed-up either with the VLA 
or with MERLIN and the VLBA to give high resolution im- 
ages (King et al., 1998). By this process the lens systems 
B0218+257 (Patnaik et al. 1993b), B1422+231 (Patnaik et 
al. 1992b) and B 1938+666 (King et al. 1997) were discov- 
ered. Also in the survey is the known system MG0414+0534 
(Hewitt et al. 1989; Hewitt et al. 1992) and the probable lens 
B2114+022 (Augusto et al. 1998). The advantage of JVAS 
compared to the MIT - Green Bank (MG) lens survey (Ben- 
nett et al. 1986) is that by observing only flat spectrum, 
hence core-dominated, sources, it is relatively straightfor- 
ward to recognize lens systems, even those which consist 
of just two images. We report here the discovery of such a 
double system, B1030+074. We present VLA, MERLIN and 
VLBA radio images and EVN information, together with 
HST WFPC2 results at 555 nm and 814 nm. In Section 2 
we present the radio data and the optical data obtained for 
B1030+074. A discussion of B1030+074 as a lens system 
follows in Section 3. 



2 RADIO AND OPTICAL OBSERVATIONS 

In this section we present all the radio maps and information 
obtained since the discovery of B1030+074. In Table 1 we 
present the radio data. Columns 1-4 show the telescope used 
for various observations, the date, the frequency and the 
resolution, respectively. The flux density in mjy of the A 
and B components as well as the flux density ratio A/B are 
shown in columns 5-7. In Table 2 we present the optical 
data. Columns 1-4 correspond to the same information as in 
Table 1 while in columns 5-7 we quote the Johnson V and 
Cousins I system magnitudes for the two components and 
their flux density ratio respectively. 



2.1 Radio observations 

All the data were analysed using a combination of the NRAO 
AIPS software package and the Caltech DIFMAP software 
package (Pearson et al. 1994; Shepherd 1997). The discovery 
map of B1030+074, made from VLA A-configuration 8.4 
GHz data giving 0.2 arcsec resolution, is shown in Fig. 1. It 
shows two distinct compact components with a separation 
of 1.56 arcsec, with the fainter image, B, at a PA of 142° 
relative to image A. The flux ratio at 8.4 GHz is 12.6. 

VLA data for B 1030+074 were also obtained at two 
more epochs, in 1994 February 22 at 8.4, 15 and 22 GHz, 
and 1995 December 20 and 19 at 15 and 22 GHz respectively. 
Though both components are detected, each is unresolved 
in the maps. The 15 and 22 GHz maps from 1995 December 
20 and 19 are shown in Fig. 1. The integrated flux densities 
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were determined using the tasks IMEAN or JMFIT/IMFIT 
within AIPS and are also shown in Table 1. 

MERLIN data were obtained at 1.7 and 5 GHz in 1993 
September 27 and 1996 December 27 respectively with reso- 
lutions of 150 milliarcsec and 50 milliarcsec respectively. The 
two components of the lens system are again unresolved and 
a model fit to the data yields a flux ratio of 18.8 at the 1.7 
GHz and 12.0 at the 5 GHz data. The MERLIN 5 GHz map 
is shown in Fig. 1. 

VLBA observations from 1995 November 12 at 5 GHz 
with a resolution of 3 mas were able to resolve at least one of 
the components. In Fig. 2 the stronger A component is seen 
to have a jet-like extension to the North-East (PA of «65° ) 
and 20 mas in length. The weaker B component remains 
unresolved even at this resolution. From the VLBA data 
the flux ratio of the two components is found to be 13.0. 

EVN observations at 1.7 with a resolution of 15 mas and 
5 GHz with a resolution of 5 mas were carried out during 
1994 May 15 and November 18 respectively. Both lens com- 
ponents are detected and are unresolved at 1.7 GHz. The 1.7 
GHz data give a flux density ratio of 18.1 for the two com- 
ponents and a position separation of 1567±1 mas (J2000) 
in position angle 143?4±0.1. However, there is a slight ex- 
tension of the A component visible in the 5 GHz map, at 
the same position angle as the VLBA extension. From the 
5 GHz data we find the flux ratio of the two components to 
be 15.9 and the position separation is the same as for the 
1.7 GHz data. 

Since we had VLA data in two epochs at each of the 
three frequencies 8.4, 15 and 22 GHz, we were able to look 
for variations in the components. There is a definite change 
in the flux density ratio and the flux densities of the two 
components change with time (see Table 1). 

2.2 HST observations 

On 1997 February 3, HST images were obtained of 
B1030+074 using the Wide Field Planetary Camera 
(WFPC2) in two filters F555W (nearly Johnson V) and 
F814W (nearly Cousins I). Two exposures of 500 seconds 
each were obtained in each filter. For the reduction of the 
data we used STSDAS and other packages within IRAF^ 
The images were rotated, cleared of cosmic ray events and 
the two exposures in each filter were averaged. The contour 
plots of the final images are shown in Fig. 3. The lowest con- 
tour level was selected to be 2a of the sky background value 
and the contour levels are separated by a factor of 2 in inten- 
sity. The HST contour plots of B 1030+074 reveal compact 
optical objects corresponding to both radio components, to- 
gether with a galaxy between them and very close to the B 
component. These data leave no doubt that B1030+074 is a 
gravitational lens system. The lensing galaxy appears not to 
have a simple smooth light distribution nor is it symmetric 
about the galaxy core. The overall extension of 1.035 arc- 
sec is in a position angle nearly perpendicular to the image 
separation; it corresponds to 6.1 kpc at the redshift of 0.599 

* IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc. (AURA) under cooperative agree- 
ment with the National Science Foundation. 



of the galaxy (see section 3), assuming Ho = 75 km sec -1 
Mpc -1 and qo = 0. Greyscale images of the two components 
and the galaxy in both the I and V filters are also presented 
in Fig. 3. The lowest intensity level is selected to be equal 
to la of the sky background value and the maximum inten- 
sity is selected such that most of the detailed structure of 
both the lens and the lensing galaxy is revealed. The transfer 
function between these two points is linear. 

As the lensing galaxy is very close to the weak B lens 
component it is not easy to separate the galaxy from the lens 
component. In order to derive accurate photometry of the 
weak component and of the lensing galaxy, we computed 
and subtracted scaled PSFs from the direct CCD images. 
By means of the Tiny Tim program (Krist 1997) we com- 
puted an oversampled numerical PSF. We then used this op- 
timal PSF in order to subtract, interactively, using programs 
within AIPS, the B component, until we saw no residuals of 
this component. More detail of the lensing galaxy is then re- 
vealed. A contour plot of the galaxy is shown in Fig. 4. The 
lowest contour level is 3a of the sky background value while 
consecutive contours differ by a factor of 2 in intensity. Al- 
though not an edge-on galaxy the lensing configuration bears 
a resemblance to that in the B 1600+433 system (Jaunsen & 
Hjorth 1997). The extended structure to the West may be 
either a spiral arm, or a second interacting galaxy. Given the 
present resolution it is difficult to say which but future ob- 
servations with the HST using the Advanced Camera might 
tell us the answer. Using the program/command "ellipse" 
within the STSDAS package in IRAF we fitted elliptical 
isophotes to the galaxy. The mean isophotal intensity as a 
function of the semi-major axis in pixels as well as the mag- 
nitude (logarithm of the isophote flux) as a function of the 
semi-major axis in pixels and the magnitude (logarithm of 
the isophote flux) as a function of the the semi-major axis 
(in pixels) to the 1/4 are also shown in Fig. 4. The profile in- 
formation seems to support the spiral galaxy interpretation 
since an t 1 ^ 4 law does not fit the surface brightness profile 
of the galaxy. Conselice (1997) show that there is a strong 
correlation between morphological asymmetry and Hubble 
type in the sense that later-type spirals show an increase in 
asymmetry. Hence, the asymmetry seen in the lens is further 
evidence for the spiral morphology. 

Photometry of the A and B components and of the lens- 
ing galaxy was performed using the AIPS task TVSTAT 
which gives image statistics within user-defined apertures. 
We then followed the standard HST photometric procedures 
(Whitmore 1997) and relied upon the values of the PHOT- 
FLAM and PHOTZPT keywords (that is the flux of a source 
with constant flux per unit wavelength in erg s _1 cm -2 A -1 , 
which produces a count rate of 1 DN per second and the 
zeropoint of the instrument respectively) appearing in the 
header of the WFPC2 frames in order to convert the flux 
densities into standard V and I magnitudes. The individ- 
ual magnitudes of the A and B components in the John- 
son V and Cousins I systems were derived from the F555W 
and F814W magnitudes. The measurements for the galaxy 
result from the integration of the flux remaining after the 
removal of the fitted PSF to the B component. The inte- 
grated I magnitude of the galaxy is 20.44+0.1. The galaxy 
is not clearly visible in the HST V image, we only see its 
compact nucleus plus some faint extended structure which 
leads to an estimate of ~22 for the total magnitude. The 
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Table 1. Radio observations of B1030+074. Absolute amplitude errors arc estimated to be 5%; the flux density ratios arc accurate to 
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Table 2. Optical observations of B1030+074. The transformed Johnson V and Cousins I magnitudes for the A and B components are 
presented and the errors are within 0.1 mag. 

Telesc. Obs. date Wavelength Rcsol. Magnitude of A Magnitude of B Flux density ratio 

(nm) (arcscc) 

HST/WFPC2 1997 02 03 555 0.045 20.34 24.10 27.4 

HST/WFPC2 1997 02 03 814 0.045 18.75 22.17 23.4 



V and I magnitudes of the A component were found to be 
20.34±0.1 and 18.75±0.07 respectively while those of the B 
component are fainter by ~3.5 mag in both filters. The op- 
tical flux ratio of the two components (A/B is 23.4 in the 
I and 27.4 in the V) is much greater than the ratio in the 
radio which is consistent with extinction due to dust from 
the lensing galaxy. The position separation of the two com- 
ponents calculated from the I HST image is 1569±2 mas in a 
position angle 143.5? ±0.1, same as that of the radio within 
the errors quoted, while the separation between the strong 
component and the center of the galaxy is 1369±4 mas in a 
position angle 142?50±0.1 and of the weak component and 
the galaxy 168±4 mas in a position angle of 156? ±0.1. 



3 B1030+074 AS A LENSED SYSTEM 

There is no doubt that B1030±074 is a gravitational lens 
system - the HST observations showing optical counterparts 
of the radio images and the lensing galaxy itself provide in- 
controvertible proof. The flux density ratio of the two com- 
ponents of B1030±074 is one of the highest for any known 
lensed system and close to the 20:1 limit adopted for the 
JVAS lens search (King et al. 1998). We attribute the differ- 
ent flux density ratios to the combined effects of variability 
and time delay. The same explanation must apply to the 
time variable flux density ratios at the same frequency. We 
plan to monitor B1030±074 to try to measure a time delay 
which we will use for a determination of the Hubble con- 
stant. 

Spectra have been obtained by Fassnacht & Cohen 
(1998) using the Low Resolution Imaging Spectrograph 
(LRIS; Oke et al. 1995) on the Keck telescope. The emis- 
sion of the background source and the lensing galaxy were 
spatially separated on the slit so spectra were extracted for 
each. From these Fassnacht & Cohen (1998) find a redshift 
for the lens of 0.599 and a redshift of the background source 



of 1.535. From the observed image splitting in the system 
they estimate the mass of the lensing galaxy within its Ein- 
stein ring to be 1.3 x 10 11 /h M© (for H = lOOh km/s/Mpc). 

We have modeled B1030±074 using a Singular Isother- 
mal Ellipsoid (SIE) mass distribution to describe the lens 
galaxy (Kormann, Schneider and Bartelmann 1994). We 
place the SIE mass distribution on the peak of the galaxy 
surface brightness distribution. Since the number of free pa- 
rameters (5 in total; PA and axial ratio of the surface density 
distribution, velocity dispersion, and source position x,y) 
is equal to the number of constraints (position x,y of im- 
age A and B, flux ratio A/B and the center of the galaxy; 
5 in total), we are able to find a mass model that repro- 
duces the image positions and their flux ratio (Table 1). The 
critical (dashed) and caustic (solid) structure of this "best 
model" are shown in Fig. 5. To investigate the stability of 
this model, we performed 10,000 Monte-Carlo simulations, 
by adding Gaussian distributed errors to the image posi- 
tions (0.3 mas), galaxy position (4 mas) and flux density 
ratio (20%) (all la errors). For each of the 10,000 models 
we solve for the mass model parameters and source position. 
We also calculate the time delay between images A and B. 
The resulting probability density distributions of the lens 
parameters and time delay are shown in Fig. 6. The error 
range on the 'best model' parameters (Table 3) indicates the 
range that contains 99% of the probability density distribu- 
tions given in Fig. 6. These ranges assume that the SIE mass 
distribution is the correct description of the lens galaxy. 

We have also tried models with an extra Singular 
Isothermal Sphere on the extension west of the lens galaxy, 
but find it has no significant influence on the position angle 
of the lensing galaxy inferred from the mass model. However, 
it does result in a smaller axis ratio and velocity dispersion 
of the main lens mass distribution. The position angle how- 
ever appears very stable. We note that the mass distribution 
is almost perpendicular to the large scale structure of the 
galaxy. 
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Figure 1. The VLA and MERLIN images of the B1030+074 system. The J2000 coordinates for B1030+074 are Right Ascension 10 h 
33 m 34.02460 s and Declination +07° 11' 26.122" . Top left: The VLA 8.4 GHz discovery map restored with a 200 X 200 mas beam. The 
contours are 0.00037 Jy per beam X (-2, 2, 4, 8, 16, 32, 64, 128, 256, 512), and the peak brightness is 0.189 Jy per beam. Top right: VLA 
15-GHz image. The contours are 0.00105 Jy per beam X (-2, -1, 1, 2, 4, 8, 16, 32, 64, 128, 256), and the peak brightness of the image 
is 0.181 Jy per beam (0.140 arcsec resolution). Bottom left: VLA 22-GHz image. The contours are 0.00144 Jy per beam X (-2, -1, 1, 2, 
4, 8, 16, 32, 64, 128, 256), and the peak brightness of the image is 0.191 Jy per beam (0.080 arcsec resolution). Bottom right: MERLIN 
5-GHz image restored with a 50 X 50 mas beam. The contours are -0.15, 0.15, 0.3, 0.6, 1.2, 2.4, 4.8, 9.6, 19.2, 38.4 and 76.8 % of the 
peak brightness value of 0.321 Jy per beam. 



The predicted time delay is around 156/^50 days with 
an error of only a few percent (assuming the validity of the 
SIE mass model). 

The mass of the lensing galaxy from the "best mass 
model" inside the critical curve is found to be 1.552 x /h 10 11 



Mq (for Ho = lOOh km/s/Mpc), consistent with the mass 
inside the Einstein radius quoted by Fassnacht & Cohen 
(1998). 
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Figure 2. The VLBA 5 GHz maps of the two components A and B of the B1030+074 system (resolution 3x3 mas). The contours are 
set to 0.15, 0.30, 0.60, 1.20, 2.40, 4.80, 9.60, 19.20, 38.40, 76.80% of the maximum value of 0.233 Jy per beam for the A component and 
2.5, 5, 10, 20, 40, 80% of the peak brightness of 0.013 Jy per beam for the B component. 



4 SUMMARY 

We have presented radio (VLA, MERLIN, EVN & VLBA) 
results and optical (HST) results for the new gravitational 
lens system B1030+074. The radio maps all show two unre- 
solved components except for the higher resolution VLBA 
map which reveals faint jet-like structure in the strongest 
component. The lensing galaxy is revealed in the HST V 
and I images. It shows substructure which indicates that is 
it not a smooth elliptical galaxy. B1030+074 is likely to be 
a lens suitable for the measurement of the Hubble constant. 



ACKNOWLEDGMENTS 

We would like to thank Dr. Tom Muxlow for reducing the 1.7 
GHz MERLIN map of B1030+074 and Dr. Sunita Nair for 
constructive comments on the paper. This research used ob- 
servations with the Hubble Space Telescope, obtained at the 
Space Telescope Science Institute, which is operated by As- 
sociated Universities for Research in Astronomy Inc. under 
NASA contract NAS5-26555. The Very Large Array is oper- 
ated by Associated Universities for Research in Astronomy 
Inc. on behalf of the National Science Foundation. MER- 
LIN is operated as a National Facility by NRAL, Univer- 
sity of Manchester, on behalf of the UK Particle Physics & 
Astronomy Research Council. We would also like to thank 
the staff of the EVN observatories and the MPIfR correlator 
for assistance with the EVN observations. This research was 
supported by European Commission, TMR Programme, Re- 
search Network Contract ERBFMRXCT96-0034 "CERES" . 



REFERENCES 

Augusto, P., et al., 1998, MNRAS, in preparation 

Bennett, C. L., Lawrence, C. R., Burke, B. F., Hewitt, J. N., 

Mahoney, J. H., 1986, ApJS, 61, 1 
Browne, I. W. A., Patnaik, A. R., Wilkinson, P. N., Wrobel, J. 

B. , 1998, MNRAS, 293, 257 
Conselice, C. J., 1997, P.A.S.P., 109, 1251 
Fassnacht, C. D., Cohen, J. G., 1998, AJ, submitted 

Hewitt, J. N., Burke, B. F., Turner, E. L., Schneider, D. P., 
Lawrence, C. R., Langston, G. I., Brody, J. P., 1989, in Lec- 
ture Notes in Physics, vol. 330, Gravitational Lenses, ed. J. 
M. Moran, J. N. Hewitt, K. Y. Lo (BerlimSpringer Verlag), 
147 

Hewitt, J. N., Turner, E. L., Lawrence, C. R., Schneider, D. P., 

Brody, J. P., 1992, AJ, 104, 968 
Jaunscn, A. O., Hjorth, J., 1997, A&A, 317, L39 
King, L. J., Browne, I. W. A., Muxlow, T. W. B., Narasimha, 

D., Patnaik, A. R. Porcas, R. W., Wilkinson, P. N., 1997, 

MNRAS, 289, 450 
King, L. J., Browne, I. W. A., Patnaik, A. R., Wilkinson, P. N., 

Marlow, D. R., 1998, MNRAS, in preparation 
Kormann, R., Schneider, P., Bartelmann, M., 1994, A&A, 284, 

285 

Krist, J., 1997, "WFPC2 ghosts, scatter and PSF field depen- 
dence", postscript document available via the STScI WWW 
page 

Oke, J. B., Cohen, J. G., Carr, M., Cromer, J., Dingizian, A., 
Harris, F. H., Labrecque, S., Lucinio, R., Schaal, W., Epps, 
H., Miller, J., 1995, PASP, 107, 375 

Patnaik, A. R., Browne, I. W. A., Wilkinson, P. N., Wrobel, J. 
M., 1992a, MNRAS, 254, 655 

Patnaik, A. R., Browne, I. W. A., Walsh, D., Chaffee, F. H., Foltz, 

C. B., 1992b, MNRAS, 259, 1 

Patnaik, A. R., 1993a, Proceedings of the 31st Liege International 
Astrophysical Colloquium "Gravitational Lenses in the Uni- 
verse", p. 311 



© 0000 RAS, MNRAS 000, 000-000 



6 E. Xanthopoulos et al. 



Figure 3. The I and V HST images of the B1030+074 system are shown in the top left and right respectively. North is up and East 
to the left. The lowest contour level is 2a of the sky background value and consecutive contours differ by a factor of 2 in intensity. The 
lensing galaxy lies closer to the weak component and appears to be extended in the almost East- West direction. We only see the nucleus 
and a little of the structure of the galaxy of B1030+074 in the V HST image. Greyscale images of the source and the lens are shown in 
bottom left and right. The greyscale minimum level is set equal to la of the sky background value. 



B1 030+074 HST I image 



B1 030+074 HST V image 
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Figure 4. A contour plot of the lensing galaxy after the faint B component has been subtracted is shown on the top left. The lowest 
contour level is 3<r of the sky background value and the contours differ by a factor of 2 in intensity. The profile of the galaxy is shown 
on the top right. The isophote mean intensity as a function of the semi-major axis in pixels appears to show a hump due to extra 
light contribution by some distorting feature in the galaxy. At the lower left and right we have plotted the magnitude (logarithm of the 
isophote flux) as a function of the semi-major axis in pixels and the the magnitude (logarithm of the isophote flux) as a function of the 
semi- major axis (in pixels) to the 1/4 respectively. 



Table 3. The "best model" and the 99% confidence intervals from the Monte Carlo simulations. 

Velocity dispersion Surface density axis ratio Position angle (N— »E) Time delay (B-A) 
(km/s) (degree) (day/h 50 ) 

258.5 (+8.4) (-5.5) 0.66 (+0.09) (-0.12) -10.2 (+14.7) (-10.1) 156.1 (+2.7) (-2.5) 
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Figure 5. The critical (dashed) and caustic (solid) structure of 
the "best" model of B 1030+074 (Table 3). The circles indicate the 
observed image positions, the cross the inferred source position. 
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Figure 6. The probability density distributions of the line-of- 
sight velocity dispersion, axis ratio and position angle of the lens- 
ing galaxy, and the time delay between lens images A and B. 
These distributions were determined by Monte-Carlo simulations, 
which include all observational errors (see text). 
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